Search for chaotic behavior in a flapping flag 

J. 0. McCaslin and P. R. Broussard 
Covenant College, Lookout Mountain, Georgia 30750 

Abstract 

We measured the correlation of the times between successive flaps of a flag for a variety of wind 
speeds and found no evidence of low dimensional chaotic behavior in the return maps of these 
times. We instead observed what is best modeled as random times determined by an exponential 
distribution. This study was done as an undergraduate experiment and illustrates the differences 
between low dimensional chaotic and possibly higher dimensional chaotic systems. 



I. INTRODUCTION 



Low dimensional chaotic behavior has been seen in many dynamical systems. The study 
of a driven pendulum^ dripping faucets,^ and many other systems have illustrated the 
transition from periodic to chaotic behavior. The dynamics of a flapping flag have also been 
claimed to be chaotic,^ but our attempt to find experimental data on flapping flags was 
not successful. There has been extensive work on the onset of flapping and flutter for paper 
and fabrics in air^ ! as well as filaments in a two-dimensional soap filmP and heavy flags in 
water However, none of these studies has characterized the dynamics of the systems studied 
as chaotic. The goal of this paper is to determine if low dimensional chaotic dynamics is 
revealed by looking at return maps of the time between successive flaps of a flag. 

II. EXPERIMENTAL SETUP AND DATA ANALYSIS 

A schematic of the experiment is shown in Fig. [T} In order to not deal with changes 
in the flag's aspect ratio due to gravity, we let the flag hang vertically and have the wind 
directed vertically. The wind was produced from a square home box fan (approximately 
0.5 m on a side) located approximately 0.75 m above the flag support rod, in order to have 
a more uniform distribution of wind. The flag was approximately 0.35 m long and 0.25 m 
wide, and the support for the flag was approximately 0.70 m above the supporting table so 
as to create a space for the air to move. Wind speed was controlled by using a variac to 
give a variable voltage to the fan. The voltage was varied from 30 to 120 V in steps of 10 V. 
The actual wind speeds were measured using a Skywatch handheld windspeed detector. Ten 
measurements of wind speed were taken for each voltage and then averaged. The relation 
between the wind speed and the voltage was approximately linear and is given in Table [T| 

The flapping of the flag was measured using a tab glued to the edge of the flag, which 
would break the beam from a HeNe laser (shown in Fig. [2]) as suggested in Ref. HI The laser's 
illumination was measured with a Vernier light sensor and Logger Pro software running on 
an inexpensive laptop. The rate of collecting illumination data was determined by making 
data collection runs at lower collection rates until the sensor began to "miss" the flapping. 
The data collection rate was set to 35 samples/s, and the data was collected for 15min for 
each wind speed. An example of the illumination versus time is given in Fig. [3j 
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To analyze the data, a prograrrPwas used to examine the illumination versus time series to 
look for the time intervals between the flag flapping, defined in our case as the time between 
drops in intensity. The code outputs a list of time intervals, denoted a.s t 1 ,t 2 , ■ ■ ■ ,t N ior each 
wind speed. These data were then used to create plots of the time intervals versus their 
order as well as return maps, where the correlation between time intervals are studied. 

III. EXPERIMENTAL RESULTS 

An example partial time series is shown in Fig. |4| The times are not uniformly distributed, 
but are more clustered around t = 0. If we plot the time intervals in a standard two- 
dimensional return map (t n+ i versus t n ) as in Fig. |5j we see that there is no apparent 
pattern in the plot, as would be expected if the dynamics controlling the flapping exhibited 
low dimensional chaos, such as seen in dripping faucet times^ or the times between swings of 
a chaotic pendulum.^ Data was collected for voltages from 30 V to 120 V (in steps of 10 V) 
(wind speeds from just under 0.5 m/s to over 4m/s). All the return maps showed similar 
behavior, as seen in Fig. |6]for the 30 V and 80 V data sets. The 30 V data exhibits a clustering 
of points not around t = 0, but instead near (t n ,t n+ i) = (1, 1). This behavior is partly due 
to the reduction of short time scales in the flapping of the flag. As the voltage decreases 
and the wind speed lowers, there are fewer and fewer rapid flapping events. In addition, 
there is a change in the distribution of times as we will discuss. The overall behavior of the 
return maps (higher density of points at low times, lower density for high times, no apparent 
structure) is unchanged. 

As mentioned, the lack of structure in the return maps is evidence for the lack of low 
dimensional chaotic dynamics in this experiment. There is a possibility that the chaotic 
behavior could be higher dimensional and not revealed in a simple two-dimensional re- 
turn map. To investigate this possibility, three-dimensional return maps were made with 
(x,y,z) = (t n ,t n +i,t n+ 2) as shown in Fig. [7] for the 120 V data set. Again, there is no ap- 
parent structure observed in this return map. We observe similar three-dimensional return 
maps for the other wind speeds. 

The return maps for voltages above 50 V (or windspeed greater than 0.8 m/s) show strong 
clustering near the origin, implying that shorter time intervals are more probable than longer 
ones. (Even for the lower values of windspeed where there seems to be a minimum time 
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between flaps, there is still clustering around the lowest times.) For the three-dimensional 
return maps, the points are also densest near the origin, and then next along the axes, with 
fewer points along the diagonals. Although there are differences in the time scaling, the lack 
of any apparent structure in the return maps motivated us to look at other possible expla- 
nations. Figure [8] shows that a histogram of the times follows an exponential distribution. 
The time constant is order 1 s, with no systematic dependence on the wind speed. For the 
lowest windspeed measured, the histogram shows a peak at times greater than zero (1.2 s 
for the 30 V data) which could be due to the natural period of the flag (estimated to be 
0.8 s). For low wind speeds the dynamics could be due to the natural periodic motion of the 
flag. However, we did not observe any indication of periodic motion in the low wind speed 
return maps. As the windspeed is increased, the peak in the histograms decreases and for 
50 V and higher the peak is located at time equal zero. 

To compare the observed histogram order, we generated a random distribution of time 
intervals with an exponential probability. A set of iV uniform random numbers Xi were 
generated in the interval (0,1), where N was equal to the number of time intervals in one 
of our data sets. Exponentially distributed numbers yi were generated from the X{ by the 
relation y(i) = — ln(xj)^ and then scaled to match the range of collected times for the data 
set considered. This set was then used to produce a simulated return map that would be 
expected if the time intervals between flaps are uncorrelated and occur with an exponential 
probability. Such a return map is shown in Fig. [9] along with an actual return map. 

By comparing the randomly generated return map to the return map for 120 V, it is 
clear that the two show very similar behavior, namely, strong concentration of points for 
small times and lower probability for large time intervals. We did not find evidence of low 
dimensional chaotic behavior for the time intervals between flaps, but instead data that 
resembles uncorrelated times with an exponential probability distribution, which would be 
equivalent to an infinite dimensional dynamic system. It is possible that we did not take 
data over a wide enough range of wind speeds. Because we did not observe any periodic 
flutter, such as seen in the heavy flag in the water experiment,^ the system should be in 
a non-linear dominated regime (even though we did see behavior at low windspeeds that 
might have been due to the natural period of the flag). We must also consider that our 
data may need to be analyzed at a higher dimensional level (looking at projections into 3D 
space of higher dimensional return maps) to verify that there is no chaotic nature to the 
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results. However, we can conclude that low dimensional chaos cannot be easily observed in 
a flapping flag experiment P 
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Tables 



TABLE I: Voltage on fan and corresponding windspeed measured at the flag position. 



Voltage ( V ) 


Windspeed (m/s) 


on i n c 

3U ± U.o 


n a n i n o o 

U.4U ± U.23 


/in i n it 

4U ± U.o 


U.52 ± U.lb 


5U ± U.o 


U.83 ± U.24 


DU ± u.o 


1 QO _|_ n 98 
1.66 ± U./o 


7U±U.5 


2.36±U.64 


8U±U.5 


2.98±U.64 


9U±U.5 


3.16±U.77 


1UU±U.5 


3.52±U.81 


110 ±0.5 


3.87±U.74 


120 ±U. 5 


4.U7±U.76 
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Figure Captions 
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FIG. 3: Example illumination change as the flag flaps. This plot is for a voltage setting of 30 V, 
or a wind speed of approximately 0.40 m/s. The plot shows two example time intervals between 
successive flaps. 
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FIG. 4: The time between successive flaps versus the sequence iteration. This data was collected 
at 120 V. 
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FIG. 5: Return map for the data collected at 120 V. 




FIG. 6: Return map for the data collected at (a) 30 V and (b) 80 V. 
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FIG. 7: Three-dimensional return map for the data collected at 120 V. The tick marks on the axis 
correspond to 1 s. 
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FIG. 8: Histogram of times between flaps for the 120 V run. 
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FIG. 9: Comparison of the return maps for times generated by an exponential probability (a) with 
the figure of data collected at 120 V as shown in Fig. |5](b). The random data was scaled to match 
the 120 V run in the number of points and range of times. 
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